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Abstract 

We embody the density matrix renormalization group (DMRG) method for 
the 19-vertex model on a square lattice in order to investigate the Berezinskii- 
Kosterlitz-Thouless transition. Elements of the transfer matrix of the 19- 
vertex model are classified in terms of the total value of arrows in one layer of 
the square lattice. By using this classification, we succeed to reduce enor- 
mously the dimension of the matrix which has to be diagonalized in the 
DMRG method . We apply our method to the 19-vertex model with the 
interaction K = 1.0866 and obtain c = 1.006(1) for the conformal anomaly. 
PACS. 05.90.+m, 02.70.-c 

I. INTRODUCTION 

The density matrix renormalization group (DMRG) method is developed to obtain eigen- 
values of Hamiltonian matrix for one- dimensional quantum systems 0. This method enables 
us to investigate a finite system with a large size by a numerical calculation of matrices, 
which can be handled within recent computer resources such as a memory size and a cpu 
speed. Recently the DMRG method has been applied to the transfer matrix of classical spin 
systems [[J . In the both cases of quantum systems and classical systems, each spin variable 
has a discrete degree of freedom. Therefore we have matrices with a finite dimension for 
these finite systems. On the other hand, the dimension of the transfer matrix becomes in- 
finite for a finite system of spins with continuous degree of freedom such as a classical XY 
model. Fortunately it is known that the XY model on a square lattice A is translated into 
a 19-vertex model for which the transfer matrix is described in terms of realization of arrow 
variables ||. An arrow variable takes three discrete states. Hence the 19-vertex model is 
called a three-state vertex model. We can construct the transfer matrix with a finite dimen- 
sion for a finite system of the 19-vertex model described by the three-state variables. The 
19-vertex model is solved by Zamolodchikov and Fateev [0] in the case that its Boltzmann 
weights satisfy the Yang-Baxter relation. This exact solution is generalized to the g-state 
vertex model by Sogo et al. 0,§J. When the Yang-Baxter relation is not satisfied, we do not 
have an exact solution yet for that system. However it is believed that a critical behavior 
of the 19-vertex model without frustration belongs to the same universality class as the 
Berezinskii-Kosterlitz-Thouless transition |||7|-|9| • 



The purpose of the present paper is to embody the DMRG method for the 19-vertex 
model and to show that the dimension of matrices, which are diagonalized in the DMRG 



method, is reduced enormously by using the ice rule of the 19-vertex model [I(J. We obtain 
a value of the conformal anomaly as c = 1.006(1) at K = 1.0866 which is regarded as the 
BKT transition point This value is consistent with a value of c expected at the BKT 
transition point. On the other hand, the value of i]/2 appears to be 0.1175(5), which is 
smaller than the expected value | at the BKT transition point. This suggests that the 
critical value of K for the 19-vertex model is smaller than 1.0866. 

In section |HJ we briefly explain the relation between the classical XY model and the 19- 
vertex model. In section |TTT| , we explain our method by which we can reduce the dimension 
of transfer matrix in the DMRG method. In section fV|, we show our results and give 
discussions about the conformal anomaly and the smallest scaling dimension. Section [V] is 
for a summary of the present study. 



II. 19-VERTEX MODEL 

The partition function of the 19-vertex model is derived from that of the XY model. 
Assigning Boltzmann weights to 19 vertices, the 19-vertex model describes the XY model 
on a square lattice for each case with frustrations or without frustrations. For spin variables, 
which take continuous values like classical XY spins or plane rotators, we have a transfer 
matrix with infinite dimension even for a finite system. Therefore the bare XY model cannot 
be handled by a numerical diagonalization of the transfer matrix. On the other hand, the 
19-vertex model is described by discrete variables which express three kinds of arrows. Hence 
we can make a transfer matrix with a finite dimension and apply the DMRG procedure for a 
finite system. In order to demonstrate the efficiency of our approach to the BKT transition, 
we apply it to the 19-vertex model without frustration for which the critical behavior of the 
usual BKT transition has to be reproduced. 

Let us briefly explain a relation between the XY model and the 19-vertex model on a 
square lattice A. The partition function Zxy of the XY model on the square lattice A is 
defined as follows: 

Zxy = II r d6k exp{K £ cos (^ - °i ~ A y)>> W 
keA J - 7T 

where i, j and k denote site indices, 9i an angle of the XY spin at site i, K an interaction 
parameter, respectively. The sum J2(i,j) is taken over all nearest-neighbor pairs of sites. 
A bond parameter related to frustration between sites i and j is expressed by Aij. The 
frustration / is defined in terms of the as follows: 

where the summation YJ P is taken over an elementary plaquette. If the value of / is a half 
odd integer, the plaquette has a frustration. 

In a region of small parameter K, we expand the exponential function in the partition 
function for the XY model given in Eq. ([[]). Hence let us start with a following partition 
function Z: 
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z = n / +7r M k n i 1 + k c ° s (^ - °j - a v)} 
= n r M k n I 1 + t ex p{^(^ - 0j - ^i)> + f ex p{- ? (^ - ^ - ^-)} I • ( 3 ) 

The integrand of the partition function Z has the U(l) symmetry as well as that of Zxy- 
The first term of the integrand in rhs of Eq.((^) is assigned to no arrow in the 19- vertex 
model. The second term is assigned to an arrow pointing from a site i to j and the third 
term to an arrow pointing from a site j to i. We use a word "arrow" even for the case of no 
arrow in a bond. Only when the value of arrows incoming to a site is equal to the value of 
arrows outgoing from the site, a weight of the arrow configuration for the site survives after 
taking integrations with respect to an angle of XY spin at the site. Otherwise the weights 
of arrow configurations do not appear in the partition function. In the six-vertex model, we 



have a similar rule to that mentioned above 10 . That rule is called as the ice rule, whose 



name comes from the property of hydrogen ions in an ice crystal. 

The 19 kinds of arrow configurations on vertices are permitted by an "ice rule", which 
is a generalization of the ice rule for the six-vertex model, when we include no arrow on a 
bond as shown in Fig.[l]. Hereafter we simply say the 19 kinds of vertices instead of the 19 
kinds of arrow configurations on vertices. The vertex weight W{vj) depends on the kind of 
the vertex Vi G {1, 2, • • ■ , 19} at site i. The value of i>j is determined by a configuration of 
four arrows as follows: 

Vi = Vi(ai,Pi,-fi,6i), (4) 

where ai, Pi, ji and Si denote arrows surrounding the site i as shown in FigfJ. Let us express 
an up and a right arrow by +1, a down and a left arrow by —1, no arrow by 0, respectively. 
For instance v j(— 1, 0, 0, +1) = 1 as shown in Fig.|l|. The ice rule is described in terms of 
ai,0i,ji and Si as follows: 

oi-i - A - 7i + k = 0- (5) 
Using the weights for 19 vertices, we can describe the partition function Z as follows: 

Z=Y,'U W (vi), (6) 

where the summation is taken over all permitted configurations of the vertex on the lattice 
A. 



III. DENSITY MATRIX RENORMALIZATION GROUP METHOD WITH 
RESTRICTIONS ON THE TOTAL VALUE OF ARROWS 

Because there exists the ice rule for the 19-vertex model, the transfer matrix, by which 
the partition function Z is expressed, becomes a block diagonal form. This block diagonal 
form is obtained in terms of classification by the value of arrows incoming to one layer 
corresponding to the transfer matrix. The layer with length L is shown in Fig|| L is the 
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width of the system. We consider the system under periodic boundary conditions without 
frustration in this and next sections. 

From the ice rule shown in Eq.(|5|), we have 



EK-A-Ti + ^O- (7) 



Because of the following relations 



ji = a i+1 , {\<i<L) (8) 
atL+i = «i, (9) 



Eq.([7p becomes 



X> = X> (io) 

1=1 8=1 

Here we note that Eq.(|9]) comes from the periodic boundary conditions. The relation given 
by Eq.(|lTJ|) means a conservation law of the total value of arrows. Let us define the total 
value of arrows N in one layer as follows: 

iv = X> (ii) 

i=i 

Now the whole transfer matrix of the 19-vertex model is classified by iV and hence has a 
block diagonal form. We will explain this in detail in the followings. 

By using the property fllQf) , we reduce an amount of calculation in the DMRG method. 
We apply the infinite system method of the DMRG framework to the 19-vertex model. 
In addition to the vertex weight W(vi) whose values are provided in Fig.[l], we define a 
renormalized weig ht W (r \vl r) ) where vf^ means a renormalized vertex and r is the number 
of a renormalization. As an initial value, we set W^^vf* 1 ) to be equal to W(vi). The transfer 
matrix for the total value of arrows, N, is composed as 

Tn\vi, fc, V3, Al|£l, <*2, £ 3 , 5 4 ) 

= E W( r \v[ r \ax,rii,a 2 ,Zi))W(v 2 (a 2 ,(3 2 ,a 3l 6 2 )) 

x W( r \vi r \a 3 , V3 ,a 4 ,£ 3 ))W(v A (a A ,fa,ax,5 A )), (12) 

(see Fig.|) where the total value of arrows N is obtained by 

N = N r ( m )+(3 2 + N r ( V3 )+(3 4 

= N r ^ 1 ) + 5 2 + N r (^)+5 4 . (13) 

We denote arrows for a vertical bond at a renormalized vertex as £j or r^. The value of arrows 
included in the renormalized vertex is denoted by N r (^i) for or N r (r]i) for rji, N r (£i) and 
N r (j]i) are equal to 5i and at the initial step of the DMRG procedure. 
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We denote an eigenvector of this transfer matrix by ippkivi, fa, Vs, AO which corresponds 

the fc-tl 
as follows: 



to the fc-th eigenvalue. As the usual DMRG method, we construct the density matrix pp k 



Notice that the density matrix is labeled by N r (rji) + /3 2 not by N. This is due to the fact 
that r/3 and (3^ are shared by two eigenvectors constructing the density matrix as shown in 
Fig.|5]. From Eq.([T3|), we obtain N r {j]i) + /3 2 = N r (£i) + 62, and therefore the density matrix 
has a block diagonal form, as well as the transfer matrix, where each block is classified by 
the total value of arrows for the half system. This property of the density matrix is the 
reason why we can reduce a dimension of the transfer matrix by our method introduced in 
the present study 

In order to construct the renormalized vertex weight W^ r+1 \ we diagonalize the density 
matrix PN r (rn)+P2,i an d obtain its eigenvectors ^v r ( 7?1 )+/3 2 ,?7i- 111 the present study, we use an 
eigenvector of the transfer matrix for the largest eigenvalue k = 1 as a target state. The 
renormalized vertex state is labeled by rf x which means that V/v r (r )1 )+/3 2 ,7y' 1 corresponds to the 
77^-th eigenvalue of PN T (rri)+fa,i- We determine the upper limit I of rf x as follows: 




(3 r+2 < m) 
> m) 



where m is the number of states taken into account for calculation of the density matrix. 
For example, at the initial step of the DMRG method, that is, r = 0, the value of I becomes 
9 which means the renormalized arrow r)[ includes two arrows. 

The last step of the DMRG method for the 19-vertex model is a construction of the 
renormalized weight for the vertex as follows: 

= EE E V Nr{m)+Ml (r] 1 ,p 2 )W ir \v{ r \a 1 ,r]i,a2,ti)) 

°2 r/i,/3 2 £i,5 2 

x W(v 2 (a2,f32,a 3 ,52))V Nr ^ l)+ 6 2 £> i (£i,5 2 ). (16) 

This last step is illustrated in Fig.|6]. The total value of arrows included in the renormalized 
vertex becomes 

N r+1 ( V [) = N r ( Vl ) + (3 2 - (17) 

We need also W <yr+1 \v^ +1 \a 3 , rj' 3 , £3)) when we return to the first step of the DMRG 
method, but we do not need to calculate it. Since in our method for the 19-vertex model 
the system has a translational symmetry, we can use W^ r+1 \vi' +1 " > ) as W( r+1 \vi r+1) ). Then 
we return to the first step represented by Eq.(|T2D in order to iterate the DMRG procedure. 
By iterating this procedure for the 19-vertex model with the restriction of the total value of 
arrows as mentioned above, we are able to make the system size L increase systematically. 
The advantage of this method is that the dimension of the transfer matrix decreases enor- 
mously by considering the conservation law of the value of arrows in this DMRG method. 
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However we need a large enough value of m in Eq.([15|) in order to obtain a good accuracy 
of numerical calculations for a large system size. The results for the m dependence of the 
present method for the BKT transition are discussed in the following section, along with 
the obtained other results. 

IV. RESULTS FOR THE CONFORMAL ANOMALY AND THE SMALLEST 

SCALING DIMENSION 

An example of the values of N r (j]i) in the case of m — 35, N = and L = 12 is shown in 
Table |. The values of rji in Table p] are put in order according to the magnitude of eigenvalues 
of the density matrix. Because the width of this system, L, is twelve, rji represents five arrows. 
Therefore the value of N r (r)i) can take one of {—5, —4, • ■ ■ , 0, ■ • • , +4, +5} as listed in the 
first column of Table [[[ Since we set m = 35, there is no eigenvalue of the density matrix 
with N r (j]i) = ±5. These are shown by crosses in Table |. When the value of arrows is —4, 
we have only one eigenvalue, which is the 31-st eigenvalue of the density matrix. This is 
expressed as iV r (31) = —4 in the third line of Table |. In the same way as for N r (r]i) = —4, 
we have two eigenvalues with N r (j]i) = —3; N r (12) = —3 and A^.(26) = —3 as shown in the 
fourth line of Table |. For other values of N r (r]i), a distribution of rji is shown in the same 
manner explained for N r (f]i) = —5,N r (f]i) = —4 and N r (r]i) = —3. The largest eigenvalue 
corresponds to iV r (l) = 0. For the case of m = 35 given in Table H, the smallest eigenvalue 
corresponds to iV r (34) = +1 and iV r (35) = —1. These two states are degenerate. In general 
those two states with N r (rji) ^ are degenerate; we have an eigenstate with N r (r]i) and an 
eigenstate with —N r (j]i) always. Therefore an eigenstate distribution of the density matrix 
is symmetric with respect to N r (r)i). 

The dimension of a matrix which has to be diagonalized in the DMRG method is now 
enormously reduced by the method mentioned in the previous section. Therefore it becomes 
possible to handle large values of m by using our computer resources. As an example we 
show reduction of the dimension of matrices of which the largest eigenvalues have to be 
calculated. In the case of N = 0,m = 35 and L = 12, the dimension of a matrix is 
3 12 = 531,441 for a simple transfer matrix method. The dimension of the matrix reduces 
down to (35 x 3) 2 = 11, 025 for the usual DMRG method The dimension of the matrix 
is now 1,545 for our method introduced in the present study. 

The largest eigenvalue of the transfer matrix belongs to the block with N = in the 
whole transfer matrix. Hence, the largest eigenvalue with N = is used to calculate the free 
energy of the system. The m dependence of the value of the free energy is shown in Fig.[7]. 
We obtain a good convergence with the value of m. In the case of L = 30, calculations with 
m = 3 14 = 4, 782, 969 give an exact diagonalization of the transfer matrix. Three kinds of 
fitting functions are used to fit the m dependence. We obtain the smallest value of standard 
deviation by using an exponential function. These results mean that the largest eigenvalue 
of the transfer matrix exponentially converges to an exact value by increasing m even at the 
BKT transition point in our method. 

The size dependence of the free energy F at K = 1.0866, which is regarded as the BKT 
transition point by Knops et al. 0, is shown in Fig.|S]. We use the following size dependence 
1 11 1 of the free energy for the system with periodic boundary conditions: 
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F/L ~ U - ^, (18) 

where denotes the free energy per site in the thermodynamic limit. We obtain the value 
of the conformal anomaly as c = 1.006(1), which is consistent with the value of c at the 
critical region. 

The second largest eigenvalue of the transfer matrix belongs to a block with N = 1. 
Hence the smallest scaling dimension which is equivalent to r]/2 is estimated by 

,/2 = = £ In (J), (19) 

where Ajj. expresses the k-ih eigenvalue in the N = n sector of the transfer matrix || . The 
critical index rj describes an algebraical decay of a correlation function in a critical region 
of interaction. At the BKT transition point, the value of rj has to be |. An interaction 

dependence of Xq is shown in Fig||. These results shown in Fig.||] are obtained by using 



m 



27 for each L. In Fig.|9], we fail to obtain the value | of Xq if we use K = 1.0866: the 
value of rj/2 obviously appears to be smaller than 1/8 at K = 1.0866. 

We investigate an m-dependence of Xq at K = 1.0866 and K = 1.0500 up to m = 120. 
Examples of extrapolations to large values of m by using the exponential function are shown 
in Fig.|H]. We define yffl by 



e^lnAffi, (20) 

Z7T 



where Ag^2 is the largest eigenvalue of the transfer matrix with the value of arrows N 
constructed by taking account of m states of the density matrix. To obtain results which 
are equivalent to an exact diagonalization, a large value of m, i.e., m = 3 19 , for the largest 
system size L = 40, which we treat. Therefore we extrapolate the obtained results yffl to 



m = oo in Fig. 10. In FigJTTj broken lines are results of fitting by using values obtained from 
m = 70, 81, 100 and 120. The dotted lines are results of fitting by using values from m = 
50, 70, 81, 100 and 120. Difference between results obtained from these two extrapolations is 
shown as errors in Fig.0. We evaluate the value of Xq by y$ — y^ and obtain the value, 
4 1} = 0.1255(1) at K = 1.0500 for L = 40. 

In Figs.|ll] and 0, we show size dependencies of Xq at K — 1.0866 and K = 1.0500, 
respectively. Every points shown in Figs.|TT] and |T2| are values estimated from extrapolations 
of m as explained above. We observe that the value of goes to finite value for large 
L. The extrapolated value is about 0.1175 and obviously less than 1/8 at K = 1.0866. 
Therefore K = 1.0866 is in the critical region. On the other hand, it is found that the size 
dependence of at K = 1.0500 is stronger than that at K = 1.0866 as shown in Fig.[T2". 
In addition, the value of Xq increases beyond 1/8 as the system size L increases. This result 
indicates that K = 1.0500 is out of the critical region. 

Thus these results at K = 1.0866 and K = 1.0500 suggest that the critical value of K of 
the 19- vertex model is in between K = 1.0500 and K = 1.0866. To determine a precise value 
of critical point K c is an interesting problem. However this problem is beyond the purpose 
of the present study. A further investigation for this point is left as a future problem. 
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V. SUMMARY 



We have succeeded to reduce enormously the dimension of matrix to be diagonalized in 
the DMRG method for the 19-vertex model by using the ice rule. It has been found that 
the m-dependence of the free energy shows an exponential convergence even near the BKT 
transition point by using our method. An accurate result is obtained for the conformal 
anomaly, i.e., c = 1.006(1) near the BKT transition point by means of the present approach. 

From the size dependence of the smallest scaling dimension x^\= r]/2), we have found 
that K = 1.0866 belongs to the critical region. However, the estimated value of r//2 is 
smaller than 1/8. On the other hand, our results for K = 1.0500 suggest that the value of 
rj/2 is larger than 1/8. These results means that the critical point of the 19-vertex model is 
in between K = 1.0500 and K = 1.0866. A further investigation for this point is left as a 
future problem. 
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TABLES 



TABLE I. An example of N r {r]i) structure in the case of m = 35, A = and L = 12. 
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FIGURES 



W(l) = K/2 W(2) = K/2 W(3) = K/2 W(4) = K/2 



W(5) = K/2 W(6) = K/2 W{7) = K/2 W(8) = K/2 



W(9) = K/2 W(10) = K/2 W(ll) = K/2 W{12) = K/2 



W(13) = K 2 /A W{U) = K 2 /A W{lh) = K 2 /A W(16) = K 2 /A 



W{\7) = K 2 /A W(18) = K 2 /4 W(19) = 1 



FIG. 1. 19 vertices permitted by the ice rule generalized to include no arrow on a bond 



FIG. 2. Four arrows surrounding a site % 



^3 — o 



1L-1 = °>L 



FIG. 3. One layer which corresponds to the transfer matrix. The value of arrows incoming to 
the layer is equal to the number of those outgoing from the layer. When a bond is shared by two 
vertices the summation for the variable on the bond is taken. 



FIG. 4. Transfer matrix composed by weights W and 
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V'S^i'^,^, AO 




V3 








FIG. 5. Construction of the density matrix by a target state V'Vfc of the system. 773 and fi\ 
are shared by two eigenvectors. Therefore we obtain a conservation rule N r (rji) + = N r (£i) +82- 



»7i 



ft 



FIG. 6. Reconstruction of VF^) in terms of eigenvectors of the density matrix with restriction 
of the value of arrows. 
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FIG. 7. m dependence of the free energy obtained by our method. We use three kinds of 
fitting functions and obtain the smallest standard deviation by an exponential function. 
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FIG. 8. Size dependence of the free energy of the 19-vertex model at K = 1.0866. Estimated 
value of the conformal anomaly from this result is c = 1.006(1). 
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FIG. 9. Interaction dependence of the smallest scaling dimension 
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FIG. 10. m dependences of at K = 1.0500 for L = 40. 
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FIG. 11. Size dependence of = rj/2 at K = 1.0866. The value of is expected to be 
0.1175(5) for large L which is obviously below 1/8 expected for the BKT transition point. The 
broken line shown in this figure is just guide for eyes. 
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FIG. 12. Size dependence of x = rj/2 at K = 1.0500. We observe a stronger size dependence 
than that at K = 1.0866 and find that the value is beyond 1/8. The broken line show in this figure 
is just guide for eyes. 
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